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Os ', Abstract 
ON 

q ■ The next-to-leading order analysis for the cross section for hadroproduction of top quark 

| pairs close to threshold is presented. Within the framework of non-relativistic QCD a sig- 

nificant enhancement compared to fixed order perturbation theory is observed which orig- 
inates from the characteristic remnant of the IS peak below production threshold of top 
quark pairs. The analysis includes all color singlet and color octet configurations of top 
quark pairs in S -wave states and, for the dominant configurations, it employs all-order 
^ ■ soft gluon resummation for the hard parton cross section. Numerical results for the Large 

Hadron Collider at - 14 TeV and *\fs = 10 TeV and also for the Tevatron are presented. 
The possibility of a top quark mass measurement from the invariant mass distribution of 
top quark pairs is discussed. 

PACS numbers: 12.38.Bx, 12.38.Cy, 14.65.Ha 



00 

o 



3 



1 Introduction 



At the CERN Large Hadron Collider (LHC) the major part of top quarks are produced in pairs. 
Due to the experience gained at the Fermilab Tevatron [1] and the huge amount of top quarks to 
be produced at LHC the reconstruction of top quarks with good accuracy will be possible [2,3]. 
A significant fraction of top quark pairs will be produced close to threshold. Thus a dedicated 
analysis of the production cross section in this region is required which is best performed within 
the framework of non-relativistic QCD (NRQCD) [4, 5]. 

The production of top anti-top quark pairs close to the kinematical threshold has received much 
attention in the context of precision measurement of top quark properties at a future Inter- 
national Linear Collider (ILC). Theoretical calculations and dedicated experimental analyses 
have demonstrated that a precise extraction of the top quark mass, its width and the strong 
coupling constant is possible [6, 7] at the ILC. The complete next-to-next-to-leading order 
(NNLO) predictions are available since many years [8]. (For earlier work see e.g [9-12].) 
Partial next-to-next-to-leading logarithmic (NNLL) [13, 14] and next-to-next-to-next-to-leading 
order (NNNLO) [15-17] predictions were evaluated more recently. 

In contrast to the linear collider, where the physical observable is the total cross section as a 
function of energy, at the hadron collider one considers the invariant mass distribution of the 
top quark pairs. Since the expected uncertainty is significant larger than the one anticipated at a 
linear collider a next-to-leading order (NLO) analysis is probably sufficient. The calculation of 
the cross section within the NRQCD framework contains as building blocks the hard production 
cross section for a top quark pair at threshold and the non-relativistic Green's function governing 
the dynamics of the would-be boundstate. Both ingredients are available in the literature since 
many years. In particular, the hard cross section for threshold tt production can be found in 
Refs. [18, 19]. In Ref. [18] the NLO formulae were derived for quark or gluon initial states 
and a quarkonium in a J pc = 0~ + color singlet state, plus possibly a parton. The general case, 
with the heavy quark system (QQ) in S -wave singlet/triplet spin state, and color singlet/octet 
configuration is given in Ref. [19], together with the corresponding results for P-waves. The 
results of Refs. [18, 19] were presented for stable boundstates. For unstable wide resonances it 
is convenient to describe the bound state dynamics through a Green's function. 

Recently a calculation of top quark threshold hadroproduction near threshold has appeared [20]. 
(For an early discussion along similar lines see [21].) The basic idea of our approach is similar 
to the one of Ref. [20]. We aim a detailed study of the top quark production based on NLO cross 
section formulae in the NRQCD framework. In our set-up all NLO sub-processes have been 
included, i.e., also those which appear for the first time in 0(a 3 s ). Furthermore, the matching 
between QCD and NRQCD as performed in Ref. [20] and the present paper is slightly different. 
Whereas in [20] the matching has been performed for the limit where the partonic center-of- 
mass energy s approaches twice the top quark mass we include the complete dependence on 
s as given in Refs. [18, 19]. Thus, formally, the result of Ref. [20] is only valid for top-quark 
production where the velocity of both quarks is small. On the other hand, in our approach the 
relative velocity has to be small whereas the top-anti-top quark system can still move with high 
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velocity. Finally, we perform a soft gluon resummation which enhances the cross section by a 
few per cent. 

Our paper is organized as follows: In the next Section details of the formalism used for the 
calculation of the NLO cross section are provided. The effects of initial-state radiation and the 
hard contribution are discussed in Section [3] and the soft gluon resummation is performed in 
Section SI The properties of the Green's function are summarized in Section [5l In Section [6] 
the building blocks are combined and numerical results for the invariant mass distribution are 
presented. Theory uncertainties due to scale variation and unknown higher order corrections 
are estimated. Summary and conclusions are presented in Section U\ 

2 The production cross section 

Let us denote the (quasi) boundstate of a top and anti-top quark with spin S and angular mo- 
mentum L by T = 2S+1 l)j' where the superscripts [1] and [8] denote the singlet and octet 
color states. The production rate is obtained from the production cross section of a top quark 
pair with invariant mass M 2 = (p t + pf) 2 and its evolution to a quasi boundstate described by 
the non-relativistic QCD. The former is a hard QCD process at a distance ~ 1 /m t and thus 
computable within the conventional perturbative expansion in a s . 

The long-distance effects responsible for the formation of a narrow boundstate are described by 
the squared wave function at the origin ^(O)! 2 or, in the language of NRQCD, by the matrix 
elements 

((X f rif,)-(^T X ))=N s N c mO)\ 2 . (1) 

Here N s = 25 + 1 and N c = 1 , 8 denote the number of spin and color degrees of freedom, respec- 
tively. We are interested in the differential distribution dcr/dM which, for narrow resonances 
with mass M n , is proportional to 6(M - M n ). For wide resonances, the case under considera- 
tion, it is convenient to convert the factor describing the sum over individual resonances into 
the non-relativistic Green's function^ 

Yw n (0)\ 2 n8(M-M n ) -> ) Im- / " v ^- = Im G(M + iT t ), (2) 
V V M n -(M + iT t ) 

with G(M + iY t ) = G^' 8 \r = 0; M + iY t ) being the Green's function at zero distance for the non- 
relativistic Schrodinger equation discussed below. Since the typical momentum scale governing 
the nonrelativistic top quark system m t v (with m t v 2 = M + iY t - 2m t , and v being the velocity of 
top and anti-top quarks) is in the perturbative regime, and the large top quark width Y t intro- 
duces an additional cutoff scale ^m t Y t , the Green's function can be evaluated perturbatively. 
As stated above the present paper is concerned with the production of top quark pairs near 

'in the case of color octet states we cannot take Eq. (0 literally but derive a corresponding formula within the 
framework of NRQCD. 
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[18] 

threshold, thus restricted to states with L = 0, i.e. T = 2S + l S j ' . The contributions to the in- 
variant mass distribution with higher angular momentum are at least suppressed by v 2 , and thus 
of higher order (beyond NLO). 

In order to obtain experimentally measurable quantities at a hadron collider the partonic differ- 
ential cross section do",y_>r/dM is convoluted with the luminosity function 



dr 



(r,/uj) = dxi f Q dx 2fi/Pi( x u^fVj/P 2 ( x 2^j)S(T-xix 2 ), (3) 

where i,j refer to partons inside the hadrons P\ and P 2 with the distribution functions fj/p { 
and fj/p 2 . The dependence on the factorization scale /Uf cancels in combination with the one 
contained in d&ij^r/dM. The differential cross section can thus be written as 

d -^L (S .M^ = V f\ T \ d ^\r^)M d ^(s,MW f ). (4) 



M — 11 (S,M Z ) = > dr 

AM ^J p L dr 



As usual s and S denote the partonic and the hadronic center-of-mass energy squared, respec- 
tively, and r = s/S . The lower limit of the r integration is given by p = M 2 /S . The partonic 
differential cross section dcr^^/dM consists of a factor F that is evaluated in perturbative 
QCD, and can be deduced from Refs. [18, 19], and a second factor, the imaginary part of the 
Green's function G'- 1 ' 8 -' 
d(X" ■ t 

M ~ Hl-^ M2 ^f) = F lHT {s,M 2 ,n 2 )—lmG [m {M + iY t ), (5) 

where the superscript of the Green's function refers to the color state of T . Eqs. © and © 
constitute our master formulae, which contain several scales and various physics contributions 
of different origin in factorized form. In particular, the soft dynamics of the parton distribution 
and real radiation is contained in the convolution of F[^j with the parton luminosity, the 
boundstate effects are described by G. Note that at NLO the Green's function G' '°^(M + iY t ) 
and the convolution of F with the parton luminosity (X ® F) are individually independent of 
the renormalization scale n r . Thus we can discuss the two parts separately in the following two 
Sections. Furthermore, it is simpler to assess the uncertainties for the individual contributions. 

Let us at this point make a comment concerning the validity of Eq. ©, which makes use of the 
NRQCD expansion assuming v «: 1 , thus being limited to the threshold region. For larger invari- 
ant masses conventional perturbation theory is applicable (see Refs. [22-24] and Refs. [25-28] 
for recent compilations of the total cross section and [29] for a proposal to measure the top-quark 
mass from the shape of dcr/dM). In the transition region the predictions from both methods are 
expected to coincide, as will be discussed below (c.f. Fig. |4]). 



3 Hard cross section 



In this Section the ingredients for the NLO corrections to the hard cross section will be collected, 
which are taken from Refs. [18, 19]. We parameterize the function Fij^j, representing the hard 
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i c [1,8] 

gg -» So 


1 C [1,8] 

gq^> s 


le [1.8] 


„„^3c [1.8] 


3 c [I- 8 ] 

gq -> Si 


3 c [I- 8 ] 

4<? -> J 5 1 


[1,5/2] 


[1,5/2] 


[3/4,6] 


[9/4,18] 


[0,32/3] 


[0, 32/3] 



Table 1: Normalization factors Njj->T for each process for [singlet, octet] color states. (N c = 3 
is used.) 



cross section for ij — » (X stands for additional partons in the inclusive cross sections), in the 
following form: 



Fjj^ T (s, M ,fi f ) = Nij^T — — 1 1 + — — C h 



3 s 



n 



x 



+ ^^(^ C (z) + ^ln C (z)) 



n 



(6) 



Here 8 ,„[i,8] = 5 - 3c [8] = 1 and zero for all other 2 — » 1 processes, and z = M 2 /S. 
quantities & c , Jl nc , and C/j all depend on i, j, and T, the functions J?l in addition on z- 



The 



The coefficients Ch originate from the hard corrections to the production process. The functions 
tt c contain the real corrections with collinear parton splitting from one of the initial partons i, j, 
and are governed by the Altarelli-Parisi splitting functions, M nc originates from non-collinear 
real emission. These individual contributions are manifest already in Ref. [18] and the appendix 
of Ref. [19], and will be listed in the following. Note, that in Eq. © we have split off the factor 
(1 + (a s /n)Ch), which we attribute to hard corrections and thus treat as a multiplicative factor 
to the terms in square brackets. 

In Tab. [T] we collect all processes of the type ij — » TX at NLO which contribute in our analysis 
and list the corresponding normalization factors Nij^r- Note that the production of a spin 
triplet color singlet state S \ via gq or qq scattering is zero up to and including NLO. This is 
because in these channels the heavy quarks are produced through gluon splitting g* — » tt, which 
is only possible if the tt is in an octet state. 

The coefficients Ch are non-vanishing only for the processes which are present also in lowest 
order [18,19]: 



Chlgg^Sf} 
Chlgg^'S™] 



C h [qq^ 3 Sf ] ] 



2 



Afi 



7T 

3 



59 
9~ 



10 ^ 16^ 
—g n f T F~^T F . 



21n2 



~4 



(7) 
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where fio = (1 1 /3)Ca - (4/3)n/7> and Cf = 4/3, Ca = 3, 7> = 1/2, ny = 5. The last term in 
Ch\qq — > 3 S [ ] ], arising from non-decoupling of the top quark in the gluon propagator, has been 
observed and discussed in Ref. [20], see also footnote 3 on page 73 of Ref. [20]. For the other 
processes hard corrections are of higher order, thus Ch is zero at NLO: 

C h [gq^ l S™] =C h [qq^ l S™] = C h [gg ^ 3 S[ hS] ] = C h [gq^ 3 Sf] = 0. (8) 



The function tt c is conveniently expressed using Altarelli-Parisi splitting functions Pijiz) intro- 
duced below [18,19] 



mgq - 

M c [qq - 

mgg - 
mgq 

Jl c [qq 



to [1,8] 
^0 



lc[l,8] 
3 



lc[1.8] 
^0 
3^ [1.8] 



3 s m 



3 S [8] 



= (l-z)PUz)\2 



ln(l-z) 
1-z 



+ 



1 



1-z 



In 



M 2 



M 2 (l-z) 2 



ZyU 



0, 

0, 



M 2 (l-zY 



C F 



a-z)p qq (z) |2 

_ ^5(1 -z) In 



z ^7 V 
ln(l-z) 
1-z 



+ r F z(i-z), 
l 



+ 



1-z 



In 



j + 



■ + C F (l-z) 



< u 2 ^ 

M 2 
\ j 



(9) 



where the conventional plus-distributior@ was employed to regularize the singularity at z = 1. 
The splitting functions Pij(z) are given by 



P qq (z) 



2C A 

C F 
T f 
2C f 



+ -+z(l-z)-2 



-1-z z 
l+(l-z) 2 



z 2 + (l-z) 2 



1 1+z 



Ll-z 



(10) 



2 The plus-distribution follows the prescription dz[ ln f(z) = Jq dz ln y\_ — \f(z) — /( 1 )], where f(z) 

is an arbitrary test function which is regular at z = 1. It is related to the p-prescription used in Ref. [19] by 
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The functions 3K m are obtained from the non-collinear contributions. For spin singlet states we 
have 



Wnc[gg^ l S [ Q 1] ] 



-C A 



12 + 1 lz 2 + 24z 3 - 21z 4 - 24z 5 + 9z 6 



6z(l-z) 2 (l+z) 3 
- 1 lz 8 + 12 (-1 + 5z 2 + 2z 3 +z 4 + 3z 6 + lz 1 ) lnz 



-C A 



12 + 23z 2 + 30z 3 - 21z 4 - 24z 5 + 9z 6 - 6z 7 



6z(l-z)(l+z) 3 L 
-23z 8 + (-12 + 60z 2 + 24z 3 + 36z 4 + 60z 6 + 24z 7 )lnz (— ) , 



-C F -(l-z)(l-lnz), 
32 C F 

32B f 
3N}. 



z(l-z), 



(11) 



where B F = (N 2 -4)/(4N c ) with JV C = 3. Note that Jl nc [gg -» ^J^] is singular at z = 1, and 
regularized by the plus -prescription. For spin triplet states one obtains 



&ncigg^ 3 S l i ] ] 



M nc [gg^ 3 S l * ] ] 



M nc [gq^ 3 S l * ] ] 



M nc [qq^ 3 S l * ] ] 



256 B F 



6C F N 2 (l-z) 2 (l+z) 3 



x 



2 + z + 2z / - 4z 4 - z 5 + 2z z (5 + 2z + z z ) lnz 



1 



108 + 153z + 400z 2 + 65z 3 - 356z 4 - 1 89z 5 



36z(l-z) 2 (l+z) 3 L 
- 152z 6 - 29z 7 + (l08z + 756z 2 + 432z 3 + 704z 4 + 260z 5 + 76z 6 ) lnz 



-Jd-z)(l + 3z) + ^i 
4 4C/r zL 



+ 2z(l +z) lnz 



(l-z)(2 + z + 2z z ) 



C F (l-z) 2 + -^-(l+z + z 2 ) 



l-z 



(12) 



The function &„ c [qq —* 3 S \^] is also defined with the plus-prescription. The leading singular 

behavior of J{ nc is given by 3\ nc {z) ~ -Ca/(1 - z)+ both for gg — » x 5q and — > J 5 j . In 
the soft limit its behavior is insensitive to the details of the boundstate and only depends on its 
color configuration. 

It is instructive to discuss the relation between the normalizations of different processes lead- 
ing to the same boundstate. For instance, the normalization A^/_>r (see Tab. Q]) for the process 
gq — » 1 Sq'^X is fixed by gg — » 1 Sq ' \ because in the collinear limit this cross section factorizes 





£®F[ij -» x 10 6 [GeV -2 ] 


£®F[ij -» T [8] ] x 10 6 [GeV -2 ] 


1 e [1,8] 
gg —* S L 

OO (J 


20.7 21.2 20.9 


63.2 62.7 60.2 


a i o 


-0.795 -1.74 -2.19 


-1.99 -4.36 -5.47 




0.00664 0.00509 0.00398 


0.0166 0.0127 0.00995 


88^S\^ 


0.175 0.127 0.0936 


6.06 4.26 3.07 






3.99 1.68 0.279 


qq^ 3 Sf ] 




23.1 23.8 23.6 


total: (^o + ^i) 


20.0 19.6 18.8 


94.3 88.1 81.8 



Table 2: The convolution £®F for LHC at the reference point M = 2m t , for the production 
of color singlet and octet states. The three columns correspond to the scale choices \x r = jj.f = 
(m t , 2m t , 4m f ). 



into the corresponding LO process and the P gq splitting function. As a consequence the cancel- 
lation of the factorization scale dependence happens among the gg and gq initiated reactions. 

Similarly, the normalization of gq — » 3 5 J is fixed by qq — > 3 S^\ In contrast, the processes 
i r i si ^ n 81 

qq — > 5 o anc * S8 ~* S i are forbidden at LO, hence the corrections have to be collinearly 
finite. In comparison to Ref. [20] the combinations Jl c + Jl n c include terms that vanish in the 
limit z — » 1. Furthermore subprocesses that appear for the first time in 0(a 3 s ) were neglected in 
Ref. [20]. The relative size of these terms will be adressed below. 

Let us now start the numerical analysis. The partonic cross sections have to be convoluted 
with the parton distribution functions (PDFs) in order to arrive at the hadronic cross section. 
We use the CTEQ6.5 [30] set for the PDFs and take a { p(M z ) = 0.118, m t = 112 A GeV and 
V5 = 14 TeV as input values. The running of a^p(ju. r ), which is the input for the partonic 
cross sections, is evaluated with the help of RunDec [31], using the four-loop approximation 
of the function. This leads to a[ 5) (ju r ) = (0.1077,0.09832,0.09050) for /u r = (m t ,2m t ,4m t ). 
Furthermore we identify renormalization and factorization scales (ju/ = fi r ). 

As stated above, the cross section factors into the convolution £<8>F and the Green's function. 
To discuss the relative importance of the various contributions individually the results for the 
subprocesses without the factor Im G(M + jT f )/m 2 are given in Tab. [21 Note that color-singlet tt 
production is dominated by by gg — » l S Q \ Color-octet production is dominated by gg — » 

plus a 25% contribution from qq — » 3 5^. The size of the remaining subprocesses (neglected 
in Ref. [20]) amounts to five to ten percent and is strongly scale dependent. The variation 
of /j. (recall /j. = fj.f = fi r ) between m t and 4m t leads to changes of £®F by +3% and +7% 
for the total singlet and octet production, respectively. In these channels the real radiation of 
partons contains large logarithmic contributions in the NLO corrections. In combination with 
the rapidly varying parton luminosity these logarithms make up for a major part of the numbers 
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quoted in Tab. [21 The origin of these large logarithms can be traced to the singular behavior 
of the cross section near z ~ 1, regularized by plus-distributions. There exists well established 
technology for the resummation of these large logarithms to all orders in perturbation theory. 
We will address this issue next. 



4 Soft gluon resummation 



The parton channels, which exhibit enhancement due to soft gluon emission are gg 



1 S^'\ and qq 



3 S (see Eqs. © and (PT21)). The relevant logarithms are contained both 
in Jl c (from initial state radiation) and ${ nc (from FSR) and read for the three leading processes: 



^thrlogVgg 
^■thrlogigg 



4C A Di-2C A ln 



M 2 



Do-f*(l 



:) In 



M 2 



^thrlogigg 



l S^]-C A D Q , 



AC F Di - 



2C F \n 



M 2 



+ C A 



Dq-^-60—z) In 



M 2 



(13) 



where Di = [ln'(l -z)/(l -z)]+ denote the plus-distributions and all ln/i 2 /M 2 parts are included 
in the definition of threshold logarithm. Whether the threshold logarithms are enhanced or not 
depends on the behavior of the parton luminosity functions near the kinematical end point r = p. 
To investigate the size of the threshold logarithms, we evaluate the contribution of the factor- 
ized hard scattering contribution convoluted with the PDFs, i.e. X<8>F separately for the three 
contributions: tree-level, singular and regular terms. (The hard corrections (1 +(a s /n)C) are 
common to all). The threshold enhanced contributions are defined in Eq. (PT31) and correspond 
exactly to the terms included in Ref. [20], while regular terms correspond to the remainder of 
&c + &lnc in Eqs. © and CGI) without plus distributions. For M = 2m t and Vs = 14TeV we 
obtain the following results 



(£®F)[gg^ l S l * ] ] 



(£®F)[qq^ 3 sf ] ] 



' 14.5 + (4.53 + 1.68)* 
14.0 + (5.66 + 1.58)* 
13.0 + (6.37 + 1.48)* 

39.3 + (16.6 + 7.26)* 
(18.8 + 6.52)* 
(20.0 + 5.83)* 

■(3.50 + 2.91)* 
(3.41+3.56)* 
(3.28 + 3.97)* 




x 10~ 6 GeV -2 , 



x 10" b GeV" 



x \Q~ b GeV" 



(14) 



The three lines correspond to fi = ft, = n r = (m t ,2m t ,4m t ). We note that in all three cases the 
contribution of the threshold enhanced terms from Eq. (TOT) is large, although the regular terms 
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in the case of qq — » 3 5 ^ are of the same order. Technically the matching applied in Ref. [20] 
corresponds to neglect all terms which vanish exactly at threshold that is for z = 1, i.e. Eqs. (fTTT) 
and (IT2]) of Section [3l The regular terms in Eq. (fl4)) . which have not been accounted for in the 
recent analysis of Ref. [20], are of the same order as the NLO sub-processes as given in Tab.[2l 

Threshold resummation proceeds conveniently in Mellin-space. To that end we calculate the 
Mellin moments with respect to z = M 2 /s according to 



F^ T (M\/u 2 f ) 



f2 ..2\ 



J dzz N - l F lHT {lM\^ f ). 




(15) 



Then, the Mellin-space expression for the threshold enhanced terms listed in Eq. (1131) read (see 
also [32,33]) 



thrlog^ 



N 

thrlog 



[qq 



3 s; 8] ] 



2C A ln / /V + C A ln/V 



Ay E -2\n 



M 2 



+C A 



2fr + 2y E 2 -2y E \n 



+ ^0 In 



= Kno^SS -> l S™] + C A lnN + C A7E , 



thrlog 

2C F \n 2 N + C F \nN 







Ay E -2\n 




V 


{»}) 



+ C A \nN 



+C, 



2{ 2 + 2y E z + -ln 









\ 




-2y £ ln 















+ C A y E , (16) 



where we have kept all dominant terms in the large-TV limit and neglected power suppressed 
terms of order 1 /N. y E is the Euler-Mascheroni constant (y E = 0.577215 . . .). 

The resummed expressions (defined in the MS-scheme) for the individual color structures of 
the hard cross sections F of Eq. © are given by a single exponential in Mellin-space (see e.g. 
Refs. [33-35]) 



F^ T (M 2 ,iuj) 
F^M 2 ^}) 



8^ T (mlM 2 f ,Mr) ^ T (mlM 2 f ,Mr) + <KN~ 1 In" N) , 



-1 1„« 



(17) 



where F\j^ T denotes the tree level term in Eq. © and the exponents are commonly expressed 
as 

In 6» HT = lnJV-gf/A) + gJ_ r (A) + . . . , (18) 

where A = /3oa s lnN/(4n). To next-to-leading logarithmic (NLL) accuracy the (universal) func- 
tions gjj as well as the functions g 2 j^ T are relevant in Eq. (fT8l . see Ref. [25] for the extension 
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to NNLL accuracy. Explicit expressions are 

A (1) 

-2- [2 - 2 ln( 1 - 2 A) + A" 1 ln( 1 - 2 A)] , 

Po 



?qq 



2 



$ Pi 



)[2A + ln(l - 2 A)] + ln 2 (l - 2,1) 



A (D /M 2 \A (1) 
■ 2-f-y £ ln(l - 2A) + In -- )-§- ln(l - 2A) + 21n 

PO \ yU, z / PO 
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where the full dependence on p r and p^. has been kept. The gluonic expressions g l gg and 8^ g ^x 

are obtained with the obvious replacement A® — » A^. The perturbative expansions of the 
anomalous dimensions are universal and well-known. We have [36] 



A ( q l) = 4C F . 



l(2) 
4 3 



A^ = 8C F 



(20) 



and all gluonic quantities are given by multiplying A® by Ca/Cf. We also give explicit results 
for the matching functions g r ^ T in Eq. (TT7T) . 



%^r[8] 



n 



2b + 2y E 2 -2y E \n 



+ 2^0 In 



1+ T c < 



2f 2 + 2y £ 2 + -ln 





-2y £ ln 





















(21) 



For phenomenological applications [37, 38] of soft-gluon resummation at the parton level one 
introduces an improved (resummed) hard cross section F res , which is obtained by an inverse 
Mellin transformation as follows, 



C+lOO 

F- r (S,M 2 ,^ } = J ^ x -^iF^ T (M 2 ,fi 2 f )-F^ T (MW f ) 



NLO 



2 ..2^ 



(22) 
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NLO 


resummed 


88 - l S™ 
qq^S™ 


20.7 
63.2 
23.1 


21.2 
62.7 
23.8 


20.9 
60.2 
23.6 


22.0 
67.8 
23.8 


23.2 
69.7 
24.0 


24.0 
70.6 
23.6 



Table 3: Comparison of the NLO and resummed result of the convolution £®F (in 10~ 6 GeV -2 ) 
for LHC at the reference point M = 2m t . The three columns correspond to the scale choices 
jj. r = jj.f = (m u 2m t , 4m t ). The NLO results can also be found in Tab. [2] 



Here is the standard fixed order cross section at NLO in QCD, while ^^Inlo 1S me 

perturbative truncation at the same order in a s obtained by employing Eq. < fT6l ). That is to say 
that for the matching we have fully expanded all formulae consistently to 0(a s ). This adds the 
hard coefficients Ch of Eq. © to the results Eqs. (fT6l) and (I2TI) . In this way, the right-hand 
side of Eq. ([221) reproduces the fixed order results and resums soft-gluon effects beyond NLO 
to NLL accuracy. 

In Section [6] we employ Eq. (1221) for phenomenological predictions by performing the inverse 
Mellin transform numerically. To that end, one should note that the treatment of the precise 
numerical matching to the exact NLO hard cross section is a matter of choice since different 
schemes lead only to differences which are formally of higher order. We have found that the 
application of the resummed result is well justified when the kinetic energy of the top-quark 
pair is a few GeV or less, see e.g. Ref. [25], where the precise numerical value is not important. 
Another issue concerns the constant terms in Eq. ([2TI) which are sometimes modified to include 
formally sub-leading (but numerically not insignificant) terms, see for instance Ref. [37,38]. 
As just explained, in the present analysis we adopt the minimal approach, i.e. we apply Eq. (f2T|) 
(including the hard coefficients Ch of Eq. ©) and account for all regular terms in Eq. CLU 
through matching to NLO. 

In Tab. [3] we compare the fixed-order NLO and resumed result of the convolution £,®F. One 
observes an enhancement up to about 10% depending on the process. 

5 Boundstate corrections 

Let us next discuss the boundstate corrections. As mentioned above, the convolution of Fjj^j 
with the parton luminosities provides the normalization of the differential cross section, while 
its shape is mainly determined by the non-relativistic Green's function. The latter describes 
the long-distance evolution of the top quark pair produced near threshold. The kinematics of 
the produced top quark pair is nonrelativistic, and the dynamics is governed by exchange of 
potential gluons leading to the formation of quasi-boundstates. The corresponding potential is 
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Figure 1: Imaginary part of the Green's functions for the color singlet (upper solid line) and 
color octet (lower solid line) cases as functions of top quark invariant mass. For comparison, 
also the expansions of G in fixed order up to 0(a s ) with (dashed) and without (dotted line) 
Y t are plotted. The imaginary part of the NNLO Green's function for the color-singlet case is 
shown as dash-dotted line. 



given at NLO by 



4na s (n r )C [h 



1 + 



An 



(23) 



with C [1] = C F = 4/3 and C [8] = C F -C A /2 = -1/6, and a\ = (31/9)C A -(20/9)7>n/. 

The color-singlet Green's function feels an attractive force, the color-octet Green's function is 
governed by repulsion and thus does not develop a boundstate. They are both defined as the 
solutions of the Schrodinger equations 

(-*'V) 2 



2m t + 



m, 



(M + iY t )}G lhS] (r;M + iY t ) = 6 {3 \r). 



(24) 



For the Green's function at zero-distance, the NLO result is known in a compact form [39] (see 
also [14]) 



G [1 ' 8] (M + iT f ) = G [1 ' 8] (r = 0;M + iY t ) = 



C^a s <4i r )m 2 t 



4n 



An 



■ gNLO + ' ' ' 
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gLO 

Snlo 



/3o\l 2 - ILL® - /c«A (1) ) + *<A (2) + ((A (0) ) 2 - 3 «A (1) - 2*^°> if/ 
+ 4 4 F 3 (i,1,1,1;2,2,1-k;i) 



L-<// 0) +Klff 



(25) 



with 



iC [1 ' 8] a,Ou r ) 



2d 



f = 



M + iT t -2m t 



(26) 



Here L = \n(ijj. r / (2m t v)) and i/^' 1 ' = 0^(1 -ac) is the n-th derivative of ^r(z) = + (J/<iz)lnr(z) 
with argument (1 -k). The Green's function in Eq. (1231 ) correctly reproduces all the NLO terms 
in NRQCD, however, it is not sufficient to describe the behavior of the Green's function in 
the vicinity of boundstate poles. It is because the exact solution to the Schrodinger equation 
has only single poles in the boundstate energy G^' ~ 1^(0)1 /(M n — M — iT t ), while Eq. (l25l) 
is an expansion around the LO boundstate poles and thus has multiple poles of a form G ~ 



l^ 0) (0)| z /(AC ; -Mf (k = 1,2 at the NLO). However, resummation of this multiple poles into 
single poles is straightforward and well-known. We refer to Ref. [39] for further details. 

In Fig. Q] we show the imaginary parts of the color singlet and color octet Green's functions in 
the threshold region. As input we use m^ s = 170.1 GeV, which to NLO accuracy corresponds 
to m, = 172.4GeV [1], and T t = 1.36 GeV [40-42]. At NLO the Green's function is separately 
renormalization scale invariant and we are free to chose fi r independent from the hard process. 
A well-motivated physical scale is fi s = m t CFa s (fi s ) = 32.21 GeV which corresponds to twice 
the inverse Bohr radius. The corresponding a s value used in Fig. |T| is a" 1 (ji s ) = 0.1401. 
It has been observed that the color-singlet Coulomb Green's function has a well-convergent 
perturbative series for this scale choice [43]. 



f (0) 



In order to see the effect of Coulomb resummation, we plot for both color states three lines: 
the full Green's function (solid line) and the expansion of G in fixed order up to 0(a s ) with 
and without top quark width (dashed/dotted). The upper three lines in Fig. \T\ correspond to 
the color singlet case and the lower three to the color octet one. The color-singlet Green's 
function shows a pronounced peak which corresponds to the tt resonance below 2m t , while for 
color octet there is no enhancement. Note that the curve for the full octet Green's function is 
very close to the one-loop expansion (taking into account the finite top quark width). Thus for 
the color octet state the Coulomb resummation effect is negligible. In addition, one more line 
(dash-dotted) for the color-singlet Green's function is plotted including the NNLO Coulomb 
potential, which is useful to estimate yet unknown boundstate corrections to the NLO color- 
singlet Green's function. As input value we again adopt the PS top quark mass [44] given above. 
Note that in the absence of full NNLO result for the Green's function and hard correction, this 
improved Green's function would not be sufficient for a full NNLO prediction. Nevertheless, the 
difference between solid and dash-dotted curves gives an indication of the intrinsic uncertainties 
of the Green's function, which is roughly 10%. 
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The expansion of G up to 0(a s ) is obtained from the gw in Eq. (|25T) as 

+ 0{a 2 s ). (27) 



-^ImG c = Im 



l H 



— — In y 

2 



In the zero-width limit (iY t — » +/0), the color-singlet curve for the expansion exhibits a step 
of height OsCf/S (for M — > 2m f ), and the color-octet curve formally becomes negative for 
v < -a s C^n/2 which corresponds to M -2m t < 0.23 GeV. Both for the singlet and octet case 
the fixed order result without Y t the imaginary part of the Green's function vanishes below 2m t . 
The qualitative difference between the solid and the short-dashed curves will be reflected in 
the comparison of our final results for the invariant mass distribution with the prediction based 
on a fixed order calculation: for the color-singlet curve we observe a sizable excess in the re- 
gion below the nominal threshold up to roughly 5 GeV above. In the color-octet case, as a 
consequence of the relative smallness of the prediction follows roughly the Born approxi- 
mation. Although the color-octet Green's function is significantly smaller than the singlet one, 
the relatively large hard scattering factor £®F for 1 Sq plus 3 S^\ which exceeds the one for 
the singlet case by roughly a factor four, quickly over-compensates the effect of the Green's 
functions. 

In the present paper we use the analytical result of the Green's function, which includes the a s 
correction (i.e. the second term in the square brackets of Eq. (1231 ) by means of the Rayleigh- 
Schrodinger perturbation approach. In Ref. [20] a numerical solution to Eq. (|24)) has been 
employed, which resums the a s corrections to all order. The numerical solution is more stable 
against scale variation and applicable over a wide range of /i,.. However, the difference between 
the two approaches is below 2% and formally of higher order. Extensive studies on higher 
order effect to the color singlet Green's function exist in the literature (see, e.g., Refs. [8,43]), 
including different implementations of the Green's function. From the experience collected in 
the linear collider studies on /f production, we expect rather large corrections from the variation 
of ji r for the color singlet Green's function of about 20% which is significantly bigger than the 
estimate from the NNLO Green's function mentioned above. In contrast to the color-singlet case 
the higher order corrections to the color octet Green's function are expected to be unimportant 
since there is no resonance enhancement and the color coefficient is small. 



6 Invariant mass distribution 

We are now in the position to combine the results of the preceeding Sections and discuss the 
cross section for the invariant top quark distribution. 

In Fig. |2]the invariant mass distributions for LHC ( VS" = 14 TeV) is shown for the three dom- 
inant processes. The bands reflect the scale variation of the convolution £,®F which for the 
color singlet case amounts to roughly ±1%. The reduction as compared to Tab. [2] and Fig. [2] 
is due to a compensation of the jx dependence after including the sub-leading NLO processes. 
Note, however, that the corresponding Green's function shows an uncertainty due to the renor- 
malization scale variation of about 20% which is well-known from top quark production studies 



15 






335 340 345 350 355 360 365 370 375 380 

M [GeV] 



Figure 2: Invariant mass distributions for leading subprocesses: gg — » S Q ' (blue and light 
green, respectively) and qq — » 3 5 ^ (green). For each process the bands take into account scale 
variation of the hard cross sections. 



in e + e~ collisions, consistent with the difference between solid and dash-dotted curves in Fig.Q] 
and thus not discussed in Fig. [2j This pattern is also evident from Fig. [3J where all production 
channels as listed in Tab.[2]are included. The width of the bands is obtained from varying renor- 
malization and factorization scales in the hard cross section as described above. The additional 
uncertainty from the Green's function, which we estimate 20% for the singlet and below 5% for 
the octet case, is not included. 

As expected, for M < 2m t the production of tt pairs is dominated by the singlet contribution. 
However, for M > 2m t one observes a strong raise of the octet contributions, in particular of 
gluon induced subprocess which for M > 2m t + 5 GeV becomes even larger than the correspond- 
ing singlet contribution. For the color-octet case the scale dependence of the hard scattering 
amounts to +7%. Considering the threshold behavior as shown in Figs. [2] and [3] it is clear, that 
the location of the threshold is entirely governed by the behavior of the color singlet (S -wave) 
contribution. Thus, as a matter of principle, determining the location of this step experimen- 
tally would allow for a top quark mass measurement, which is conceptually very different from 
the one based on the reconstruction of a (colored) single quark in the decay chain t —> Wb. In 
fact, much of the detailed investigations of tt threshold production at a linear collider were per- 
formed for this particular relations between the location of the color singlet quasi-boundstate 



16 



3.5 I — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — r 






335 340 345 350 355 360 365 370 375 380 



M [GeV] 

Figure 3: Invariant mass distribution including all production channels shown in Tab. [21 The 
width of the bands reflect the scale dependence of the hard scattering parts. 



pole of tt and the top quark MS-mass. The absolute normalization of the cross section is also 
sensitive towards electroweak corrections [45^1-9] which are of the order of 5% close to thresh- 
old. For example, the difference between corrections from a light (Mf, = 120 GeV) and a heavy 
(M h = 1000 GeV) Higgs boson amounts to roughly 6% [48]. 

In Fig. |4] the prediction for dcr/dM based on NRQCD is compared with the one obtained 
from a fixed order NLO calculation for stable top quarks which is obtained using the program 
HVQMNR [50]. As expected from the comparison of solid and dotted curves in Fig.[T] the two 
predictions overlap for invariant masses around 355 GeV. Above 355 GeV relativistic correc- 
tions start to become important. From this comparison we find an additional contribution to 
the total cross section for //production of roughly 10 pb, which could become of relevance for 
precision measurements. Note that the band of the NRQCD-based prediction only contains the 
uncertainty from the scale variation of £,®F whereas the one of the Green's function (which 
can reach up to 20%, see Section [5]) is not shown. 

The analysis of this work has concentrated on the threshold region and is applicable for M up 
360 GeV at most. However, it is obvious, that the overall shape of dcr/dM will be distorted and 
the mean (M) shifted to smaller values, which might affect the global fit of dcr/dM. In Fig. |5] 
we present for comparison the NLO prediction for dcr/dM in the wide range up to 700 GeV. 
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Figure 4: Invariant mass distribution dcr/dM from NRQCD and for a fixed NLO for LHC with 
yfs = 14 TeV. The bands are due to scale variation from m t to 4m t . For the NRQCD prediction 
the additional uncertainty due to the Green's function estimated to 20% (5%) for the colour 
singlet (octet) contribution is not included. 



The distribution reaches quickly its maximum of 3.3 pb/GeV at around 390 GeV and then falls 
off slowly. It is remarkable that its value at 370 GeV is already not too far from the maximum 
of the curve and the threshold modifications thus affect a sizeable part of the distribution. 

Although the most detailed top quark studies will be performed at the LHC at an energy of 
14 TeV, a sample of top quarks has been collected at the Tevatron in proton anti-proton colli- 
sions at 1.96 TeV. Furthermore the first LHC data set will be taken at 10 TeV. For this reason 
we give the results for these two cases, in Figs.|6]and|71 The cross section in Fig.[6]has the same 
characteristic shape as the one in Fig. [31 however, the absolute size is considerably smaller. As 
expected, the enhancement at threshold is significantly less pronounced for Tevatron where the 
colour singlet contribution is very small. 

Our analysis confirms the findings of Ref. [20], however, the numerical results for the cross 
sections as presented in Fig. [3] are slightly higher than the corresponding corrections of Ref. [20] 
which is due to the combined effect of the soft-gluon resummation, the inclusion of the NLO 
sub-processes and the different matching to full QCD. 
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Figure 5: Invariant mass distribution dcr/dM from NLO calculation for LHC with V* = 14 TeV. 

7 Summary 

A NLO analysis of top quark production near threshold at hadron colliders has been performed. 
The large width of the top quark in combination with the large contribution from gluon fusion 
into a (loose bound) color singlet tt system leads to a sizable cross section for masses of the 
tt system significantly below the nominal threshold. A precise measurement of the M t j dis- 
tribution in this region which is dominated by the color singlet configuration could lead to a 
top-quark mass determination which does not involve the systematic uncertainties inherent in 
the determination of the mass of a single (colour triplet) quark. Furthermore, also the shape 
of the differential distribution dcr/dM is distorted and the mean (M) shifted towards smaller 
values. 

The effects of initial state radiation as well as boundstate corrections are taken into account in 
consistent manner at NLO. As compared to Ref. [20] we include the complete s dependence 
in the matching condition and also implement all NLO sub-processes. We observe a partial 
numerical cancellation between these two effects leading to similar predictions as Ref. [20]. 
Furthermore we perform a soft-gluon resummation and thus include the dominant logarithmi- 
cally enhanced higher order terms. This last step stabilizes the prediction. However, it enhances 
the cross section at most by 10%. 
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Figure 6: Invariant mass distribution dcr/dM for LHC with V* = lOTeV. 



The effects are more pronounced at the LHC with top production being dominated by gluon 
fusion and less relevant in proton-antiproton collisions with top quarks dominantly in color octet 
states. Considering the threshold region (say up to M t j = 350 GeV) seperately, an integrated 
cross section of 15 pb is obtained, which should be compared to 5 pb as derived from the NLO 
predictions using a stable top quark and neglegting the binding correction. Within this relatively 
narrow region the enhancement amounts to roughly a factor three and a significant shift of the 
threshold. Compared to the total cross section for //production of about 840 pb (obtained using 
fixed-order NLO accuracy for n = m t , see, e.g., Ref. [25]), the increase is relatively small, about 
1%. However, in view of the anticipated experimental precision of better than 10% these effects 
should not be ignored. 
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Figure 7: Invariant mass distribution Act I AM for Tevatron with ^[s = 1.96 TeV. At the Teva- 
tron qq — » 1 Sq dominates the cross section, and luminosity for gg channels is small, thus the 
boundstate peak is buried by color-octet production. 



Note added 

While this article was finished an analytic evaluation of the total cross section at NLO accuracy 
appeared [51], which has been used in Ref. [20] to clarify the existence of a non-decoupling top 
quark effect overlooked in Ref. [19] (see footnote 3 on page 73 in Ref. [20]). 
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